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Some features of the Driven Jaynes-Cummings system 

B. TugulduiQand Ts. Gantsog 

School of Physics and Electronics, National University of Mongolia 

We investigate a generalized Jaynes-Cummings model with the external field driving the cav- 
ity mode. The numerical results for the dynamics of the atom and the cavity field mode are 
given. Approximated solutions to the Schrodinger equation are supplied, which allow us to 
predict the kinematics of the Q(a, a*) function in the complex plane. 



1 Introduction 

The Jaynes-Cummings model (JCM) HI is one of the simplest models describing the interaction 
of light with matter, where a single two-level atom interacts with a single mode of quantized radi- 
ation field in the electric dipole and rotating wave approximations. This model is of fundamental 
importance to the field of quantum optics; many interesting features have been predicted for both 
the atomic variables and the statistical properties of the field through the years, beginning with 
the well-known phenomena of collapses and revivals of the atomic population inversion oscilla- 
tion 0. Numerous extensions of the JCM have been considered OHH and many experiments 
reported ll6tiT0l. 

One of the most interesting features of the JCM is that if the atom is initially prepared in its 
upper state and if the cavity field is initially in the coherent state, then the quantized field evolves 
into an almost pure state at half of the atomic-revival time lfTTHT3l . This approximately pure 
state is equal to a superposition (Schrodinger cat) state composed of two states of light having 
the same amplitude, but opposite phase. The amplitude of the component states is approximately 
equal to the amplitude of the initial coherent state of the field mode lfl4l . 

Buzek et al. lfT31l showed that by driving the atom with the external classical field, superpo- 
sition states of the quantized cavity mode with arbitrary amplitudes and phases of component 
states can be produced. Recently, Gea-Banacloche and coworkers |fT6l have verified the exis- 
tence of maximally entangled state in externally driven JCM. Optical cavities with atoms have 
been proposed for quantum information processing 111711181 . 

In this paper we study the dynamics of the JCM when a quantized cavity mode is pumped 
continuously by an external classical field. In Sec. 2 we derive the eigenstates and eigenenergies 
of the driven Jaynes-Cummings system. Using these results we find the time dependent state 
vector of the system for a given initial condition. In Sec. 3 the dynamics of the atomic inversion, 
mean photon number and the phase space distribution function Q(a, a*) are examined. In Sec. 4 
we give approximate solutions that enable us to explain the features observed. Finally in Sec. 5 
we summarize our results. 
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2 Driven Jaynes-Cummings System 



We consider the driven Jaynes-Cummings system with the external field driving the cavity mode. 
The interaction Hamiltonian of the system in the interaction picture is given by [3] 

Hi = g(aa + + tf&J) + <?(aV + fie - **), (1) 

where g is the coupling constant between the atom and the cavity mode; S is the amplitude of 
the driving field; <fr is the phase of the classical field; a', a are creation and annihilation operators 
for the cavity mode; and <x + and <r_ are atomic pseudospin operators. We assumed that the 
driving field is in resonance with the atom and the cavity mode. The steady state solution of the 
Schrodinger equation for the Hamiltonian (1) is possible for § < g/2. The quasieigenenergies 
and the corresponding eigenstates are given by J5] 

£o = 0, 



E n i = g\JTi 
E„2 = -g\fn 
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and 



\4>Eni) = c p (S( V )D(l3 n+ )\n - 1) + L+S( V )D(p n+ )\n))\-) 

+ c p (L + S( V )D(f3 n+ )\n - 1) + e 2i *S(r))D((3 n+ )\n})\+), (3a) 

|V> B „ 2 ) = c p {S(ji)D(fi n J)\n - 1) - L + S{r,)D{p n _)\n))\-) 

+ c p (L+S(r,)D(p n _)\n - 1) - e 2 ^S(?7),D(/3„_)|n))|+), (3b) 

where if quasienergies are positive (negative) E n \ > (E r a < 0) then j3 n+ = \fne 1 ^ 
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= ^\-{2Sjgf, 



= 2(f), 



(1/2)^/1 - y/1 - {2S/g) 2 and L+ = -{ge^/2S)(l + yjl - {2S/g) 2 ). For an external field 
amplitude larger than <§ = g/2, no normalizable steady states exist. This value is the threshold 
condition for spontaneous dressed-state polarization. 

If the atom is initially prepared in its upper state |+) and if the cavity field is prepared initially 
in the coherent state \ao), then the time evolution of the system is given by 

\m) = w + (t))+\i>-(t)), (4) 

where 



H + (t)) = K\ 2 X>+> - !; -f?|a ) + e-*+{n; -p n+ ;- v \a ))e 



iE n t 



[(L+\ m (3 n+ ;n - 1) + e 2l %- /?„+; n))\+) + (\ m (3 n+] n - 1) + L+\m /3„+; n))|-)] , 

(5a) 
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OO 

n=0 

[(L+h; /?„_; n - 1) - e 2l *\ m /?„_; + (\r,; /?„_; n - 1) - L + \ m /?„_; n))|-)] , 

(5b) 

and £?„ = <7v^(l — (2<£7<?) 2 ) 3 ^ 4 and squeezed and displaced Fock states are denoted as 
S(f])D((3 n ±)\n) = 1 77 ; j3 n ±;n). This result is exact and no approximations are made. We shall 
apply it to study quantum dynamics of the system in the following section. 

3 Quantum dynamics of the driven JCM 

3.1 Atomic inversion and the mean photon number of the cavity mode 

Using explicit expressions for the state vector given by equations and © we can study the 
dynamical properties of the system under consideration. Firstly we calculate the expectation 
value of the atomic inversion. An illustration of the time evolution of the atomic inversion os- 
cillation in the driven JCM is shown in figure [TJ for different values of § . The revival of the 
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Fig. 1. The atomic inversion as a function of gt for the atom initially in the upper state and the cavity field 
in the coherent state |«o). The parameters are ao = 5, g = 1, <j> = and (a) £ = 0.01; (b) S = 0.05; (c) 
S = 0.1; (d) S = 2. 

atomic inversion oscillations is a purely quantum-mechanical effect that originates in the dis- 
creteness of quantum states of the cavity mode. It is seen from figure [TJ that with the change 
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of the amplitude S of the driving field the revival time of the atomic inversion as well as the 
overall pattern of the time evolution of the atomic inversion oscillation are changed. For larger S 
the overall revival and collapse pattern disappear due to the domination of the classical driving 
field on the quantized cavity field. In the standard JCM {§ = 0) the total excitation number 
operator R = a^a + 1/2(1 + a z ) commutes with the total Hamiltonian and hence represents 
an integral of motion. As a consequence, the time evolution of the mean photon number of the 
cavity mode shows the same oscillation as the atomic inversion. However, in the case of driven 
JCM ($ ^ 0) we have [R, H] =^ 0, and therefore the total excitation number is not conserved. 
Physically this corresponds to the fact that the energy is transferred from the classical field to the 
quantized cavity mode. Using equations (01 and ([5]) we can numerically calculate the mean pho- 
ton number of the cavity mode. In figure [2] we plot the mean photon number (a^a) for ao = 5, 
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Fig. 2. The mean photon number as a function of gt for the atom initially in the upper state and the cavity 
field in the coherent state \ao). The parameters are ao = 5, g = 1, <j> = and (a) $ — 0.01; (b) S = 0.05; 
(c) S = 0.1; (d) S = 2. 

= and for different values of S '. As we can see from figure |2]the time evolution of the mean 
photon number in the driven JCM exhibits several interesting features. Firstly, for small values 
of § (see figure [2(a)| with S — 0.01) a typical collapse-revival pattern corresponding to the time 
evolution of the inversion can be observed. At the same time we can also observe the pattern of 
the global oscillating behavior. With the increase in the classical amplitude S the time evolution 
of the mean photon number exhibits stronger global oscillating behavior with small quantum os- 
cillations (revivals) at times corresponding to revivals of the atomic inversion (the amplitude of 
these oscillations is of the order of unity and therefore they are not clearly seen in our pictures 
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Fig. 3. The atomic inversion as a function of gt for the atom initially in the upper state and the cavity field 
in the coherent state |qo). The parameters are ao = 5, g = 1, $ = 0.1 and (a) <j> — (when <j> — ir, the 
result is exactly the same as this one); (b) 4> = n/2 (the vertical bars indicate the times in Fig. 5). 



because of the scale used in the y axis). Finishing this part of the paper we turn our attention to 
the figure [3] where the atomic inversion is plotted for «o = 5, & = 0.1 and for different values 
of 0. We can see from the pictures that the revival time of the atomic inversion as well as the 
overall pattern of its time evolution also depend on the phase of the driving classical field <j>. 

3.2 Q-function 

Now we turn our attention to the dynamics of the Q(a, a* ) function. We know from the standard 
JCM (S = 0) that, if the initial field is in a coherent state \cto), then during the evolution of the 
system the single-peaked Q-function of the cavity mode bifurcates into two peaks that move in 
opposite directions around a circular path whose radius equals the amplitude ao of the initial 
coherent state ifTTl . The centre of this circular path is at the origin of the phase space. At one 
half of the revival time the two peaks have a maximum separation in phase space. Phoenix 
and Knight lfl2l and Gea-Banocloche fT31 have shown that at this time the cavity mode (and 
the atom) returns most closely to the pure state, which is close to a macroscopic superposition 
composed of two field states that have the same amplitude but opposite phase. The purity of this 
macroscopic superposition state increases as the amplitude of the cavity mode increases. As the 
peaks meet at the opposite side of the circle, they produce a revival of the Rabi oscillations. 

Buzek et al. ifTSIl have shown that, by driving the atom with the external classical field, one 
can produce (almost pure) macroscopic superposition states of the cavity field. The amplitudes 
of the component states can be much larger than the amplitude of the initial cavity field mode. 
With a proper choice of the phase of the classical driving field, one can produce superposition 
states with different amplitudes and phases of the component states. 

It is interesting to contrast this behavior with what occurs when the cavity quantum mode is 
driven by an external classical field. In figure |4]the contour plots of the Q(a, a*) function are 
shown for ao = 5, S jg = 0.1 and for three different values of the phase <f> = 0, ir and ir/2. The 
most interesting result is that in the course of the time the initial one-peaked function splits into 
two peaks, which counterrotate on separate circles with different radii depending on the phase 
of the classical driving field <f>. It is clearly seen from the comparison of figure [3(b)1 and figure[5] 
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that atomic inversion oscillation shows revival when these two peaks collide. When = and 
<f> = it, the two circles have different radii, and the peaks can collide only at the original site of 
the circle (see figures |4(a)| and |4(b)| . But when S /g is small enough, the peaks still do have some 
overlap while moving on different circles (not necessarily on the opposite sides of the circles) 
and as a result we observe some small revival pattern (see the first small revival in figure [Kb)) . 
However, with the increase of S /g, the peaks do not overlap and the revival pattern disappears 
(see figure [T(c)] >. 

In the next section we will derive an approximate expression for the state vector that enables 
us to study the kinematics of the Q(a, a*) function in the complex plane: we will be able to 
determine radii and centers of two separate circles, angular velocities of rotation of the peaks 
around these circles, and therefore estimate the revival time of the atomic inversion oscillation. 



4 Approximate solutions 

If <§ is small compared to <?, we can approximate the state vector dD-© of the system to give 
(see Appendix A) 

\i,+ (t)) « A+(t)\(ao + Pn)e- lu>+at - n )(\+) + |-», (6a) 



|VT(t)> « A- (t)|(a -A,)e & '""* + &)(|- 
where f3 n = {2S'/g)Vne ief ', and 



(6b) 



A+(t) = | Cp | 2 |L + | 2 exp(-|ao+^| 2 



exp \ -l-P„(a + Pn)*e w gt + ^(a + ^ n )e^ 9t ] \ exp 



(7a) 



A-(t) = |c p | 2 |L + | 2 cx P ( -|a -/3„|' 



cx P \ -(/3„(a - /9n)*e-* w gt - fc{a - p n )e w ~*) } exp 



and (see Appendix B) 

u+ [1 - (2^/g) 2 ] 3 / 4 
2\a +p n \ 

u - = [1 - (2^/g) 2 ] 3/4 
2|ao-A,| 



(7b) 



(8a) 



(8b) 
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Fig. 4. Contour lines of Q(a, a* ,t) in the complex a plane for ao = 5, g = 1, § = 0.1 in the initial 
condition - |qo)|+), at the times gto — 0, gt\ = 11.5, gt-2 — 23, gt$ = 34.5. (a) (j> = 0; (b) cj> — tt; (c) 
(f> = vr/2. 



It is clear from the above results that ?/» + corresponds to the coherent state rotating along the 
circle of radius |arj + f3 n \ centered at — /3 n , and %p~ corresponds to the coherent state rotating 
along the circle of radius |ao — /3 n | centered at f3 n . These states counterrotate in the complex 
plane. This approximation is valid only if time is less than revival time. 
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Fig. 5. Q(a, a* , t) function of the cavity mode for the times gto = 0, gti — 21, gti = 42, gt-j, = 51, 
gti — 67. The parameters are the same as in Fig. 3b. The times are indicated by the vertical bars in Fig. 3b. 



The modules of the A + and A are not time dependent. Therefore density operator of the 
cavity mode takes the form 

Pf = Tr atom {p) = 2|A+| 2 | 7 +(i))( 7 +(<)| + 2| J 4-| 2 | 7 -(t))( 7 -(t)|, (9) 

where p is the density operator for the total system, 7 + (t) = (ao + j3 n )e~ luJ+9t — /3„, and 
7 _ (£) = (ao — fin) e%u 9t + Pn- Then Q(a,a*) function can be written as a sum of two 
counterrotating peaks in the phase space i.e., 

Q(a,a*,t)^|H7 + W)| 2 + |(«|7"W>| 2 - (10) 

In figure [6] we plotted the time evolution of Q(a,a*,t) using exact, as well as approximated 
solutions for = (figures [6(a)| and |6(c)| i and for (f> = n/2 (figures [6(b)| and |6(d)| i. The initial 
condition of the system is j+) jao). The left hand side pictures correspond to the exact solutions 
and the right hand side pictures correspond to the approximated results. The initial Q(a,a*) 
distribution is a Gaussian centered at ao. As we know, in the course of time the initial one- 
peaked function splits into two peaks. We plotted them in separate Figures. One of them (V> + ) 
rotates counterclockwise (see figure [6(a)) on the circle of radius \cto + 1| = 6 (/3 n ~ 1) with 
angular speed uj + = 0.0808, and the other one (i/^ - ) rotates clockwise (see figure |6(c)] i on the 
circle with smaller radius |ao — 1| = 4 at angular speed u~ = 0.1212. In this case the centers of 
circles are shifted respectively down and up in real axis. For (f> = ir/2, ifj + and counterrotate 
on the different circles with the same radius |arj±i| = 5.099 at angular speed cj ± = 0.0951. The 
centers of the circles are shifted horizontally. We can see from figure |6]that our approximated 
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Fig. 6. Contour lines of Q(a,a* ,t) in the complex a plane for the initial condition |+)|ao) and for 
ceo = 5, g = 1, $ = 0.1. (a) (j> — 0, ip + (t); (b) </> = 0, tp~ (t); (c) and (d) are the approximated results and 
parameters are same as (a) and (b) respectively. 
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Fig. 7. Contour lines of Q(a,a* ,t) in the complex a plane for the initial condition |+)|ao) and for 
ceo = 5, g = 1, $ = 0.1. (a) = 0, ip + (t); (b) </> = 0, tp~ (t); (c) and (d) are the approximated results and 
parameters are same as (a) and (b) respectively. 
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solutions give almost precise locations of the peaks. Therefore we can predict the kinematics of 
the moving peaks and find the times when those peaks collide, which enable us to calculate the 
revival times of the atomic inversion. 



5 Conclusion 



We have calculated steady states of the generalized Jaynes-Cummings system with the external 
field driving the cavity mode. The steady states are the superposition of atomic states multiplied 
by squeezed and displaced Fock states. The calculated steady states are valid only if 2S/g < 1. 
In other words, 2S /g = 1 is a threshold value of the driving field strength and above this value 
normalizable steady states do not exist. But we numerically calculated the time evolution of 
normalized states for the arbitrary value of 2S '/g. Using these results we studied the dynamics 
of the atomic variables and the statistical properties of the field. 

Furthermore, we derived approximate solutions for the time evolution of quantum states, 
which allow us to understand the kinematics of Q(a, a*) function in the phase space. We found 
that as time goes on the initial single peaked Q(a, a*) function splits into two peaks which 
counterrotate on separate circles of different radii \ao ± /3„| with different angular velocities u> ± . 
Whenever these peaks collide, or even partly overlap, atomic inversion oscillation shows revival. 



A Approximation method 

If we are interested only in the kinematics of the Q(a, a*) function in the complex plane, but 
not the shapes of the peaks, then we can make some very simple approximations for the state 
vectors (0. We restrict our calculations with only first order terms of S '/g ignoring all higher 
order terms to yield r « 0, and therefore rj rj 0, and S(r)) rj 1 (squeezing operator does not 
affect the location of the peaks in phase space, it changes only the shape of the peaks). 

Doing same approximation we get L + w — (g/S)e 1 ^. Since \L+\ is much greater than l,the 
terms not containing L + can be neglected in equations ©. Then we get 

00 

H' + (t)) « K\ 2 \L + \ 2 ^«n-l;/? n |ao»e- iE " 4 (| - p n \n- 1)|+) + | - /?„;n>|-» , (11) 



\r{t)) « M 2 |L+| 2 £«n - l;-p n \a ))e iE ^ (\p n] n-l)\+) - |0„;n>|-» , (12) 

n=0 

where f3 n = (2<? / g) y/ne 1 ^ . The difficulty to manage above equations is that (3 n depends on the 
photon number n. The multiplier IS / g in equation (3 n is much less than 1. That leads to (3 n ^ 
in small interval of n. If we change j3 n with 

[3 n = (2S/g)Vne^, (13) 

we can easily evaluate state vector 

\^{t)) « |c p | 2 |L + | 2 £ {ao+ J. n)n e-^D(-p n ) (\n - 1)|+) + \n)\-)) , (14) 

n=0 
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\^-(t)) « |c p | 2 |L + | 2 f; K e^D{p n ){\n-l)\+) \n)\-)) . (15) 



If ao is large enough, then «o + Pn and ao — (3 n are also large and there is not much difference 
between \n — 1) and \n). We can substitute e Tlul nt instead of e^ %Ent . Therefore, taking into 
account the relation 



D(a)D(/3) = cxp 



D(a + /3), (16) 



we get the approximate expressions for the state vectors as 



\^(t)) « A+(t) £ [(a ° +/ ^" )e ^ +gt hT \n) (|+) + |-» , (17) 

n— 



|r (*)> « £ ^^^ In) (|+) - I-)) . (18) 

n— 

B Estimations of w+ and cj 

Here we make a simple estimation for the rotational angular velocities oj + and oj~ of two separate 
peaks of the Q(a,a*) function in the complex plane. First we estimate oj + . Using equation ( fT4b , 
the reduced density operator for the quantized cavity field can be written as 

p+ g (a + &,)" ((a e _, fcgt(VH _^ g( _^ )|w)H d[ - Mj (19) 

n.rn— 

where fc = (1 - (2<f /g) 2 ) 3 / 4 . Using a Taylor expansion of ^/n around \/n (where n = |ao + 

V 7 ?! « =(n — fi) — (20) 

and ignoring all higher order terms we get 

sfri — \frri m [n — m)/2A/n. (21) 
Substituting it into the equation ( fT9l ) for one can find w + as follows 

, , [1 - (2<f / 5 ) 2 ] 3 / 4 



2|a +£„| 

Similarly, for we get 

[1 - (2^/g) 2 ] 3 / 4 
2|a -^n| 



(22) 



(23) 
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